Metal-insulator transition of the Kagome lattice fermions at 1/3 filling 
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We discuss the metal-insulator transition of the spinless fermion model on the Kagome lattice at 1/3-filling. 
The system is analyzed using exact diagonalization, the density-matrix renormalization group methods and the 
random phase approximation. In the strong-coupling region, the charge-ordered ground state is consistent with 
the predictions of the effective model with a plaquette order. We find that the qualitative properties of the metal- 
insulator transition are totally different depending on the sign of the hopping integrals, reflecting the difference 
of band structure at the Fermi level. 
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Introduction — For a long time, the exotic phenomena 
emergent in geometrical frustration have been fascinating 
but challenging subjects of research in condensed matter 
physics iit]. An essential difficulty of frustrated systems is 
attributed to their highly-degenerate ground states in the clas- 
sical sense. It is common practice to determine how the 
degeneracy is removed or how the frustration is minimized 
by adding quantum fluctuations. On a related issue, one of 
the most intensively studied systems is quantum spins on 
the Kagome lattice, which consists of corner-sharing trian- 
gles. Very recently, so-called herbertsmithite has been re- 
vealed as a promising candidate for an experimental real- 
ization of an S" = 1/2 Kagome antifeiTomagnet. Both ex- 
perimentally 10, [3] and theoretically i4[|, the low-temperature 
phase of this system is characterized by a spin liquid state with 
a small or possibly even no energy gap. 

When we turn our attention to itinerant systems on the 
Kagome lattice, a greater variety of physical phenomena can 
be observed. Of particular interest is the competition of dis- 
tinct ordered states associated with the charge degrees of free- 
dom. From this point of view, an intriguing system is that 
of 1/3-filled Kagome lattice fermions. If long-range repulsive 
interactions are taken into account, there exists a huge num- 
ber of (nearly) degenerate charge-ordered (CO) states in the 
strong-coupling regime; a melting of the CO may take place 
as some ratio of the interaction and hopping. 

An experimental candidate to realize this situation is 
hydrogen-bonded crystals of alkali-hydrosulfates or hydrose- 
lenates M3H(X04)2 (M=Cs,Rb,K,Tl, X=S,Se) H. These 
materials exhibit an abrupt change in the electrical conduc- 
tivity accompanied with a ferroelastic deformation around 
Tc = 400K: above the conductivity within the a-b plane 
is very high and it seems to be related to the hopping motion 
of protons, i.e., super-proton conductivity; below T^., the pro- 
tons appear to be localized in some ordering pattern which 
satisfies the condition that one tetrahedron of XO4 has only 
one hydrogen bond (see Fig. [T]). In the regime where only 
proton ordering is taken into account, an effective model can 
provide a description of the protons forming the Kagome net- 



work, as discussed later in the paper Since the number of 
corner-sharing triangles is equivalent to the number of tetra- 
hedra, the critical behavior in these materials can be described 
by a metal-insulator transition (MIT) in the Kagome lattice at 
1/3 fining. 

Another possible experimental realization may be given 
by recently developed laser-cooling techniques. It was re- 
ported 101 that a two-dimensional optical trinierized Kagome 
lattice was constructed and that the coupling constants could 
be controlled using a triple laser beam design. Furthermore, a 
scheme proposed this year would construct an optical Kagome 
lattice at 1/3 filling with the use of just two standing waves and 
more controllable interactions ^ . 

Motivated by the optical lattice systems, a model of hard- 
core bosons on the Kagome lattice was studied with large 
scale quantum Monte Carlo simulations and phenomenologi- 
cal dual vortex theory At 1/3 and 2/3 fillings a weak 

first-order superfluid-solid transition induced by the nearest- 
neighbor repulsion was found. Thus, in our work presented 
here it is natural to consider how this transition occurs in 
fermion systems. Although the effective model can be shown 
to be equivalent for hard core bosons and spinless fermions 
in the strong-coupling limit, features for each, in the weak to 
intermediate coupling regimes, may differ completely due to 
the absence of a band picture in boson systems; here the dis- 
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FIG. 1 : (Color online) Lattice structure of hydrogen-bonded crystals 
M3H(X04)2 where the proton sites foiTn the Kagome lattice (right 
figure). The left figure presents an example of protons order below 
the critical temperature. The solid (red) lines denote hydrogen bonds. 
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FIG. 2: (Color online) Band structure of the Kagome lattice fermions 
for positive (left) and negative (right) hopping energies. 



tinction between the two cases is non-trivial. 

In this Letter, we therefore study the MIT of the spin- 
less fermion model on the Kagome lattice at 1/3 filling 
with the nearest-neighbor repulsive interactions using the 
exact-diagonalization (ED) and density-matrix renormaliza- 
tion group (DMRG) techniques. Assuming a CO state in the 
strong-coupling regime, we demonstrate that the nature of the 
MIT strongly depends on the sign of hopping integral (inver- 
sion of the sign corresponds to a switch between 1/3 and 2/3 
fillings). Lastly, the random phase approximation (RPA) is 
employed to confirm the numerical results. 

Model — The Hamiltonian is 



{id) 



(1) 



where c] (cp is a creation (annihilation) operator of a spinless 

fermion and jij {— c^jCj) is the corresponding number opera- 
tor. The repulsive interaction V{> 0) is assumed to act only 
between neighboring sites In the non-interacting case 

V =^ 0, the dispersion relation is given by three bands: 

e{k) = -2t, 
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(see Fig. It is particularly worth noting that the flat band 
at 1/3 filling is only filled when t > 0. Since the density of 
states at the Fermi level diverges, the physical properties in the 
weak- to intermediate-coupling regime are expected to greatly 
differ from those for t < 0. We note that our model with t > 
at 1 /3 filling is equivalent to that with t < at 2/3 filling via 
the particle-hole transformation. 

Let us now consider the case of t = with V > 0. The 
nearest-neighbor repulsions V are minimized when each the 
corner-sharing triangle is occupied by exactly one fermion. 
This means that the system is in macroscopically degenerate 
CO state. However, the degeneracy may be lifted when a finite 
hopping is introduced, and then a CO state with some period- 
icity must have the lowest energy. The candidates of the CO 
patterns are quoted in Fig. [3] By looking at the pattern III, 
one realizes the kinetic energy can be gained through a cyclic 
hopping process of three fermions on each hexagon. Thus, 
the lifting takes place with processes of order |<p/y^ ( lower 



order contributions, e.g. t'^ /V, are the same for all configu- 
rations and therefore contribute a constant energy shift rather 
than lifting the degeneracy). The effective Hamiltonian is de- 
rived as uBUll 
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(3) 



where the sum is over all hexagons on the lattice and the six 
sites on each hexagon are labeled by 1-6 in a clockwise man- 
ner This is equivalent to the quantum dimer model on the 
honeycomb lattice which has a plaquette-ordered ground state 
with three-fold degeneracy 11311 . The effective model would 
be at least qualitatively valid as long as the system is a CO in- 
sulator We note that here the sign of hopping integral has no 
influence on the physics because the honeycomb lattice has a 
bipartite structure. 

Numerical analysis — First, we compare the lowest ener- 
gies of each CO pattern to determine if pattern III is indeed 
stabilized for y ^ |t| > 0. We apply the ED technique with 
periodic boundary conditions. A 12- (18-)site cluster is used 
to calculate an energy difference Ei — Eu {En ~ Em) [see 
Fig. m, where Ei, Eu, and i^iii are the energies of the CO 
patterns I, II, and III, respectively. The periodicity of wave 
functions is checked in order to identify one CO state from 
another. As shown in Fig. [5] for both signs of t, Em has 
the lowest ground-state energy while £"11 is the second-lowest 
in the large-F regime; the energy differences are scaled like 
Eu - Em 0{V~^) and Ei ~ Eu 0(y-^). Since 
the second-order contributions Odtj/y^) cancel in the en- 
ergy differences, the behavior of £"11 — £"111 seems to be con- 
sistent with the existence of an energy gain (!?(|i|'^/y^) from 
the ring exchange process in the effective Hamiltonian (O. As 
expected, £11 — £111 has no dependence on the sign of t in the 
large-y region as seen in Fig.|5](b). 

Next, we show that the system indeed exhibits the CO- 
melting MIT at a finite critical interaction V = Vc- To find 
the MIT critical point, the single-particle gap is calculated us- 
ing the DMRG method. The DMRG calculation is performed 
on a Lx X Ly cluster shown in Fig. lUb), where the periodic 
boundary conditions are imposed in the y direction with three 
unit cells and the open boundary conditions in the x direction 




FIG. 3: (Color online) Charge-ordering patterns in the 1/3-filled 
Kagome lattice characterized by wave vectors I {qa.,qh) = (0,0), 
II (0, tt), and III {2n/3, 2tt/3). The circle on pattern III corresponds 
to the cyclic hopping process by Eq.ljS]!. 
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with Lx unit cells. We study several lengths of the cluster with 
up to La- = 8 (78 sites) and extrapolate the finite-size results 
to the thermodynamic limit L^. oo. The single-particle 
gap is thus given by A = \vai]^^^^ ^{L^) with /S.{Lx) = 
E{Nf + 1, i„ Ly) + E{Nf - iX, Ly) - 2E{Nf,Lx, Ly), 
where E{Nf,Lx,Ly) is the ground-state energy of the cor- 
responding cluster with Nf fermions. The upper panels of 
Figure |6] show the extrapolated values of A as a function of V 
in units of \t\. The critical points are estimated as Vc ^ 2.6 
and 4.0 for t > and t < 0, respectively. The behaviors 
of the gap opening are quite different: the gap for <> in- 
creases almost linearly, while it rises gradually for t < like 
the Berezinskii-Kosterlitz-Thouless transition. 

In order to gain further insight, we also calculate the or- 
der parameter corresponding to the CO pattern III. The CO 
phenomenon is observed as a state with a broken transla- 
tional symmetry: actually, there are two or more degenerate 
ground states and one configuration of the degenerate states 
is picked out as the ground state by an initial condition of 
the DMRG calculation, as a consequence of the cylindri- 
cal boundary conditions [14J. The order parameter is now 
defined as the amplitude of the charge-density modulation 
rj{Lx) = Til — \{i^2 + "-a) for a x 3 cluster, where ni 
and 712 (= "^3) are the charge densities of fermion-rich and 
fermion-poor sites on a triangle at the center of the cluster. 
In the lower panels of Figure |6] we show the extrapolated re- 
sults of order parameter 77 to the thermodynamic limit, i.e., 
rj = limi^^oo il{Lx), as a function of V . We find that the 
estimated critical points are almost the same as those obtained 
from A for both signs of t. As seen similarly for A., the be- 
havior of 77 differs markedly depending on the sign of t. For 
t > 0, the transition appears to be nearly discontinuous, which 
might be analogous to that in a melting MIT observed in hard- 
core bosons on the triangular lattice 11511 : for t < 0, it looks 
continuous, which looks like a weakly first-order phase tran- 
sition for hard-core bosons on the Kagome lattice ||8|]. Hence, 
we conclude that the criticality for positive t is in stark con- 
trast to that for negative t. 

Now, we comment on the difference in values of the MIT 
critical points between i > and t < 0. As mentioned above, 
all the triangles are only allowed to be singly occupied in the 
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FIG. 4: (Color online) (a) Finite clusters of the Kagome lattice ap- 
plied for the ED method. The 12- and 18-site clusters are used to ob- 
tain El — Ell and £11 — -Em, respectively, (b) La: x 3 cluster (shaded 
regime) for the DMRG method, where the open and periodic bound- 
ary conditions are taken for the x and y directions, respectively. 
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FIG. 5: (Color online) (a) Energy differences £"11 — Em as a function 
of in units of \t\, and (c) Ei - En as a function of (c) 
Relationship of the energies Ei, En, and Em corresponding to the 
CO patterns I - III in Fig. [3] 




FIG. 6: (Color online) Single-particle gap A and order parameter 77 
as a function of the nearest-neighbor repulsion V for (a) positive and 
(b) negative t values. Insets: Finite-size scaling of each quantity for 
(a) V = 2, 3, 4, and 5; and (b) V = 2,2,, 4, 5, and 6, from bottom to 
top. 



CO state. Therefore, our system can be mapped into a single- 
band half-filled Hubbard model on honeycomb lattice with 
on-site Coulomb interaction Oiy) if we regard each trian- 
gle as a single site. The single-particle gap is then expected 
to behave like A w I^ — VFcff in the large-y region, where 
Wcff is the effective bandwidth. Actually, as seen in Fig. |6] 
the gap behaves like A - F - 5t and A - ^ - 8|t| for t > 
and t < 0, respectively. In general, the MIT critical strength 
is roughly scaled by the bandwidth, i.e., I4 Woff. The as- 
sumption of the effective bandwidths Wcff = ht for i > and 
%\t\ for t < are compatible with the obtained critical points. 
The narrower bandwidth for t > might be interpreted as a 
remnant of the flat band in the noninteracting case ||T6[| . 

Random phase approximation — In addition, the MIT is in- 
vestigated analytically using the RPA iH. The charge sus- 
ceptibility is given by the Kubo formula as 

1 t^'' 

X^^{q,iOJl) ^ — dTe'"'^(71q,,(T)71_q,^(0)) (4) 

-'''uc JO 
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FIG. 7: (Color online) Phase diagram of the 1/3-filled t-V model on 
the Kagome lattice obtained by the random phase approximation. 



where ( • • ■ ) is the ensemble average, riq „ is a number opera- 
tor in momentum space, /3 = /cbT is inverse temperature, 
is the Matsubara frequency of bosons, and N^c is the number 
of unit cells. Then, applying the RPA, the charge susceptibil- 
ity is obtained in the following matrix form, 

Xiq,iLJi)^ [l + 2z;(q)X(o)(q,ic^0)]"'^^°^(g,i^/), (5) 
where the matrix elements of v{q) are given by 

Vi2{q) = Vcos{qr,/2), Vi3{q) ^ V cos{qy/2) (6) 
V23{q) = Vcos{{q,,-q.y)/2), (7) 



with the relation v 
ity is defined as 



3 1 



and Vii = 0. The bare susceptibil- 
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where 



ihuJl — £k+q,v + Efc,/ 



(9) 



with / being the Fermi distribution function. Uk is a matrix 
which diagonalizes the kinetic part of eq. ([T]l in the momentum 
space k. 

The MIT critical point at finite temperature is found by de- 
termining the divergence of the diagonal part of eq. Q for the 
three sets of CO patterns, I, II, and III. We found that the CO 
pattern III is always most stable, in agreement with the ED 
result. The obtained phase diagram is shown in Fig. [7] At 
high temperature {k'QT/\t\ ^ 1), the MIT phase boundary is 
almost independent of the sign of t; however, for decreasing 
temperature the boundaries gradually move apart from each 
other The critical points are estimated as 14 ~ 3 and 4.5 for 
t > and t < 0, respectively. These values seem to be not 
only qualitatively but also quantitatively consistent with those 



obtained by the DMRG method. It is also interesting that the 
boundaries show reentrant behaviors. 

Conclusion — We study the metal-insulator transition of the 
spinless fermions on the Kagome lattice at 1/3-filling using the 
exact diagonalization, density-matrix renormalization group, 
and random phase approximation techniques. In the region of 
large V, the CO pattern in the ground-state is consistent with 
what is predicted by the ring exchange model. The behaviors 
of the single-particle gap opening for positive and negative 
hopping integrals are qualitatively quite different; this is fur- 
ther reflectedin the differing values of the critical point Vc- 
These differences may arise from the distinct nature of the 
density of states at the Fermi level in each case. A possible 
further extension of this work is a consideration of spin de- 
grees of freedom. This is of interest because spins in the CO 
phase are also frustrated and the configuration is nontrivial. 
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